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1.  INTRODUCTION 

The  first  work  on  the  dynamical  theory  of  the 
oscillation  of  liquid  drops  appears  to  have  been 
done  by  Lord  Rayleigh.'  In  his  investigation, 
Rayleigh  considered  small  displacements  from 
equilibrium  under  the  action  of  forces  due  to 
surface  tension  (capillary  forces)  alone.  Later, 
Lamb2  included  damping  of  the  small  oscillation  of 
drops  by  the  inclusion  of  internal  viscous  forces 
and  showed  that  the  rate  of  damping  was  depen¬ 
dent  on  the  size  of  the  liquid  drops,  becoming 
extremely  large  for  very  small  drops.  The  possi¬ 
bility  of  instability  of  the  lower  modes  of  oscilla¬ 
tion  was  first  investigated  by  Rayleigh,3  who 
showed  that  the  resonance  frequency  of  the  nth 
mode  of  a  charged  drop  was 


n(n  —  1) 


y(n  +  2) 


where  p  is  the  density  of  the  liquid  drop,  a  is  the 
equilibrium  radius,  y  is  the  surface  tension,  and 
Q  is  the  total  charge  on  the  drop.  Instability 
occurs  when  wj,  <  0;  or  when  the  charge 
]Q  |  >  [4na3y(n  +  2)]l/2,  then  the  nth  mode  is 
unstable;  the  lower  modes  become  unstable  for  the 
smaller  amount  of  charge.  The  condition  for  insta¬ 
bility  is  sometimes  written  in  terms  of  the  electric 
field  at  the  surface  of  the  drop,  E  =  Q/a2,  so  that 
the  condition  expressed  above  is  frequently  given 
as  |E|  >  [47r(n  +  2)y/a]1/2  and  has  been  used  as 
a  starting  point  by  later  investigators. 

Considerably  later  than  the  above  work,  the 
investigation  of  the  dynamics  of  small  droplets  was 
begun  by  Thacher4  and  O’Konski  and  Thacher.5  In 
their  investigations,  they  considered  the  distortion 
of  droplets  by  an  electric  field,  assuming  that  the 
distorted  drop  was  an  ellipsoid  of  revolution  about 
the  direction  of  the  electric  field.  They  also  dis¬ 
cussed  the  possibility  of  enhancing  the  droplet’s 
distortion  by  an  alternating  electric  field,  but  failed 
to  observe  that  the  frequency  of  the  electric  field 

*Lord  Rayleigh.  Proc.  Roy.  Soc  (London),  29  (1879).  71.  The  Theory  of 
Sound.  II.  MacMillan  Co..  London  ( 1896).  Ch.  XX. 

*H.  Lamb.  Hydrodynamics.  Dover  Publications.  Inc..  He *■  York  (1932) 

3 Lord  Rayleigh.  Phil.  Mag..  14  (1882).  184. 

*H.  Thacher.  Jr..  J.  Phys.  Chem..  56  (I9S2).  79 5 

5C.  T  O’Konski  and  H  Thacher.  Jr  .  J  Phys  Chem  .  51(1953).  955 


should  be  half  that  required  for  resonance  of  the 
drop.  O’Konski  and  Harris6  extended  this  work  to 
include  in  their  analysis  the  effect  of  electrical 
conductivity  of  both  the  surrounding  medium  and 
the  drop.  No  consideration  was  given  to  the 
possible  effect  of  conductivity  on  the  damping.  In 
their  analysis,  they  found  the  rather  surprising 
result  that,  under  certain  appropriate  choices  of 
conductivities,  the  equilibrium  shape  of  the  droplet 
remained  a  sphere.  As  in  their  previous  investiga¬ 
tion,  O’Konski  and  Harris  assumed  that  the  drop¬ 
let  under  the  action  of  an  electric  field  becomes  an 
ellipsoid  of  revolution  about  the  field  direction. 

Later,  Garton  and  Kras;ucki7  investigated  the 
stability  of  bubbles  in  a  static  electric  field  and 
questioned  the  correctness  of  previous  theoretical 
work  on  electrostriction.  They  showed  actual 
photographs  of  bubbles  in  various  stages  of  dis¬ 
integration  that  vividly  displayed  the  physical 
nature  of  bubbles  breaking. 

In  a  series  of  two  papers,  Taylor8  discussed  the 
stability  of  conducting  drops  in  an  electric  field.  He 
showed  that  it  is  necessary  to  introduce  motion  of 
the  fluid  inside  a  drop  to  attain  the  spherical 
(undistorted)  solution  of  O’Konski  and  Harris.6 
Some  experimental  work  was  reported  in  Taylor’s 
papers  that  agreed  quite  well  with  the  theory. 
Sozou9  later  extended  Taylor’s  theory  to  include 
time  dependent  electric  fields.  He  gave  a  number  of 
results  in  the  form  of  equations  with  a  few  cases  of 
numerical  results,  but  further  work  seems  neces¬ 
sary  for  comparison  with  experiment  and  with 
previous  results. 

Rosenkilde10  investigated  the  stability  of  drop¬ 
lets  in  an  electric  field  using  methods  of  tensor 
calculus  and  Chandrasekhar’s"  virial  method.  He 
showed  that  under  appropriate  conditions  at  least 
three  different  equilibrium  configurations  can  exist. 
He  predicted  instability  only  if  the  dielectric  con¬ 
stant,  r,  of  the  drop  is  greater  than  20.801. 

T  O  'Konskt  and  F  E  Harris.  J.  Phys.  Chem..  61  (195  7),  1172. 

7G  G  Gorton  and  2.  K.asucki,  Proc.  Roy.  Soc.  (London),  A280 (1964).  211. 

8 Sir  Gtoffm  To- tor.  Pro  \  Roy.  Soc.  (London).  KUO  ( 1964),  383;  A291 
(19*6).  159 

V  Swim.  Proc  Roy  Soc  (London).  A331  (1972).  263. 

E  /(osenktlde.  Proc  Roy.  Soc.  (London).  Kill  (1969).  473. 

*  *S  <  nandrasekhar.  Hydrodynamics  and  Hydro  magnetic  Stability.  Claren¬ 
don  Prr  .s,  Oxford  (19*1) 
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Most  of  the  above  work  was  devoted  to  the 
development  of  equilibrium  distortion  of  a  droplet 
under  the  action  of  various  forces.  None  of  these 
workers  attempted  to  develop  their  results  using 
the  original  techniques  introduced  by  Rayleigh. 1  In 
fact,  the  technique  used  by  Rayleigh  in  his  dis¬ 
cussion  of  a  charged  drop3  was  somewhat  different 
than  in  his  original  paper.  A  detailed  derivation  of 
Rayleigh’s  result  for  a  charged  drop  was  given  by 
Hendricks  and  Schneider12  and  later  by  Schneider13 
along  with  many  applications  of  the  results.  Dis¬ 
sipation  by  viscosity  was  not  considered  in  either 
of  these  investigations. 

In  our  consideration  here,  we  follow  the  spirit 
of  the  original  investigation  by  Rayleigh.  In  section 
2. 1 ,  we  find  the  energy  due  to  the  surface  tension, 
the  kinetic  energy  of  the  fluid  in  the  drop,  and  the 
Rayleigh  dissipation  function.  From  these  quan¬ 
tities,  we  derive  the  generalized  equations  of 
motion  by  using  the  Lagrangian,  as  given  by 
Rayleigh,  which  includes  the  results  of  Lamb. 

In  section  2.2,  we  include  the  terms  involving 
the  energy  due  to  the  electric  field  in  the 
Lagrangian.  We  discuss  this  aspect  in  consider¬ 
able  detail  since  we  have  not  been  able  to  find  any 
reference  to  an  extension  of  Rayleigh’s  result  that 
includes  an  electric  field. 

The  remainder  of  the  report  is  devoted  to  a 
discussion  of  the  result  of  section  2.2.  The  results 
are  used  in  the  investigation  of  instability  under  the 
action  of  an  electrostatic  field  as  well  as  the 
dynamical  behavior  of  drops  under  the  action  of  a 
time  varying  field. 

2.  THEORY 

2. 1  Derivation  of  Rayleigh’s  Result 

In  the  dynamical  theory  of  the  oscillation 
of  a  liquid  drop,  we  assume  that  the  radius,  r,  from 

^  Lord  Rayleigh,  Proc  Rov.  Soc.  (London).  29  (1879).  71;  The  Theory  of 
Sound.  !L  MacMillan  Co..  London  (1896),  Ch.  XX. 

^Lord  Rayleigh.  Phil.  Mag..  14  (1882).  184. 

,2C.  D.  Hendricks  and  J.  M.  Schneider.  Am.  J.  Phys..  31  (1963).  450. 

M John  Matthew  Schneider.  The  Stability  of  Electrified  Liquid  Surfaces. 
Charged  Particle  Research  Laboratory.  Vnivet  ty  of  Illinois.  Urbana.  CPRL-2-64 
(5  March  1964 )  (NT1S  604431). 


the  center  of  mass  to  a  point  on  the  surface  of  the 
drop  can  be  expanded  in  a  Legendre  series  as 

r(0,t)  =  ao(t)  +  £'  ak(t)  Pk(cos  0)  .  (1) 

k 

The  method  then  of  solving  any  particular  problem 
is  to  express  the  Lagrangian  for  that  problem  in 
terms  of  the  variables  a„(t)  and  treat  these  a„(t) 
as  generalized  coordinates  to  obtain  the  equations 
of  motion.  Such  an  expansion  as  equation  (1)  is 
possible  if  we  assume  that  the  droplet  is  symmetric 
about  the  z-axis.  At  this  point,  the  z-axis  can  be 
chosen  in  any  direction;  but  later,  when  we  include 
the  electromagnetic  energy,  the  electric  field  will 
be  assumed  to  be  along  this  axis.  We  assume  this 
symmetry  throughout  the  discussion.  The  prime  on 
the  sum  in  equation  ( 1 )  indicates  that  the  lowest 
value  of  k  in  the  sum  is  one. 

With  the  assumption  of  incompressibility, 
the  volume  of  the  drop  acts  as  a  constraint  on  the 
ak(t)  in  equation  (1).  The  volume  of  the  drop  is 
given  by 

X2jt  pi i  pi(fiX) 

d </>Jb  sinfldfljf,  x2dx  .(2) 

Since  we  are  assuming  that  the  drop  is  symmetric 
about  the  z-axis,  letting  y  =  cos  0,  we  have 

27T  Z'1 

V  =  — r3(M)  dy  .  (3) 

Using  equation  ( 1 )  and  the  relation 

S-\  P*(M)P*”(M)  dy  =  28gg'/(2g  +  1)  , 

we  obtain 

V“  y(«  +  3-?2irfr)  •  <4> 

If  we  assume  that  a  is  the  radius  of  the  equilibrium 
sphere  (V  =  4n-a3/3),  then 

■»”•('  'I?'?'  27Tt)  • 
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which  holds  for  terms  through  order  ak  with 
k  >  1.  Since  we  are  interested  only  in  terms  of 
second  order  in  ak,  we  treat  equation  (5)  and 
corresponding  subsequent  results  as  equalities. 

The  potential  energy,  Us,  of  the  drop 
due  to  the  surface  tension,  y,  is  the  area  of  the 
drop  multiplied  by  y  or 

us  =  y  Jo  ^ S  rsin^‘ls  .  (6) 

where  ds  is  the  arc  length  in  the  surface  given  by 


V  •  V  =  0  .  (11) 

If  further  we  assume  that 

V  X  v  =  0  ,  (12) 

then  7  can  be  derived  from  a  potential  function, 
0, ,  such  that 

v  =  -V0i  ,  (13) 

and  from  equation  ( 1 1 )  we  obtain 


ds2  =  dr2  +  r2  dO2 


V20,  =  0  .  (14) 


or 


ds  =  dO 


Thus,  the  velocity  potential  satisfies  Laplace’s 
equation,  and  the  solution  appropriate  for  our 
(7)  problem  can  be  written 


Using  this  result  in  equation  (6),  we  can  write 


0,  =  £  /J„rnP„(cos  6)  .  (15) 

n 


Us  =  2ny 


/.fI 


r2  +  (1 


(8) 


The  second  term  in  the  square  root  in  equation  (8) 
is  smaller  than  r,  and  the  square  root  can  be 
expanded  and  integrated  to  give 


Us  =  4  iry  I  a 


a 


+  r 

k 


(k  -  lXk  +  2)4 
2(2k  +  1) 


L  >(9) 


where  the  constraint  given  in  equation  (5)  has  been 
used  to  obtain  this  result.  The  result  given  in 
equation  (9)  holds  up  to  fourth  order  in  the  ak. 


To  calculate  the  kinetic  energy,  T,  we 
need  to  evaluate  the  volume  integral 

T  =  -y  y*pv2  dr  ,  (10) 

with  p  the  density  (constant)  and  v  the  velocity 
of  the  fluid  within  the  drop.  We  assume  that  the 
fluid  is  incompressible  and  that  there  are  no 
sources  or  sinks;  then 


(Here  r  represents  a  point  on  the  interior  of  the 
drop  and  becomes  r(0,t)  at  the  surface.)  Rather 
than  use  equation  (15)  in  equation  (10)  directly, 
we  use  equation  (13)  in  equation  (10)  to  convert 
the  volume  integral  into  a  surface  integral  as 

T  =  yp  /0,(V0.)  ■  do  .(16) 

We  assume  that  the  area  element  do  is  approxi¬ 
mately  along  r.  (Actually,  do  =  (rrdfl  —  0  dr) 
X  r  sin  9  d 0,  where  r  and  0  are  unit  vectors 
along  r  and  6,  but  the  correction  is  of  higher 
order  than  we  are  considering  here.)  Under  that 
assumption,  we  have 


ka2k~' 

T  =  2opa2  S’  yyy  Pi  .  (17) 

The  Pk  in  equations  (15)  and  ( 1 7)  can  be  evalu¬ 
ated  by  equating  the  velocity  at  the  surface 
(—d<pjd r)  with  r  calculated  from  equation  (1) 
so  that 

ak  =  -ka  k_,/?k  (18) 

Then  equation  (17)  becomes 


7 


”  ■ T-  MB 


T 


2npai 


£' 

k 


»k 


k(2k+  1) 


(19) 


4npa3 
n(2n  +  1) 


an  +  8; rr] 


a(n  —  1)  . 

— : - &n 


This  result,  equation  (19),  is  valid  only  through 
terms  of  order  a£.  If  the  result  given  in  equation 
(18)  is  carefully  inspected  in  the  relation 

Vr  =  -  £  k/^-'P^cos  0) 


4ny(n  —  l)(n  +  2) 
2n  +  1 


an 


=  0 


(22) 


Thus,  if  the  viscosity  vanishes  (rj  =  0),  we  have 
Rayleigh’s  result,1 


and  the  rk-1  is  expanded  by  using  equation  (1), 
correction  terms  occur  in  fa.  These  terms  give 
corrections  to  in  the  form  a^a,,,  thereby 
coupling  the  equations  of  motion  of  in  a  com¬ 
plicated  manner,  which  we  ignore  in  the  present 
investigation.  However,  any  serious  investigation 
of  dynamic  instability  should  include  these  terms 
because  they  are  of  third  order  in  the  and  can 
be  important. 


w n(b  =  0)  =  -b-  n(n  -  l)(n  +  2)  .  (22a) 
pa 


If  the  viscosity  is  small  enough,  we  obtain  Lamb’s 
result2  for  the  decay  time  of  the  nth  mode 
[an  an(0)  e_t/rnj. 


pa2 

r,(n-l)(2n+l) 


(22b) 


To  derive  the  effect  of  viscosity  on  the 
equations  of  motion,  we  use  Rayleigh’s  dissipative 
function,14  R,  instead  of  the  procedure  used  by 
Lamb.2  As  will  become  apparent  in  the  develop¬ 
ment,  this  method  is  much  more  consistent  with  the 
treatment  presented  here.  A  convenient  form  of  the 
dissipative  function  given  by  Landau  and 
Lifshitz13  is 

R  =  y  r,  f(S7  v2)  •  da  ,  (20) 

where  rj  is  the  viscosity.  As  in  the  kinetic  energy, 
it  is  sufficient  to  assume  that  da*  is  along  7.  All  of 
the  quantities  necessary  to  calculate  equation  (20) 
are  given  in  equations  (15)  and  (18)  so  that 

R  =  4Trrja  £'  ~  ~  .  (21) 

k  K 

The  equation  of  motion  for  a„(t)  can  now 
be  readily  found  by  forming  the  Lagrangian  from 
the  results  given  in  equations  (9)  and  (19) 
(L,  =  T  —  Us)  and  the  dissipative  function 
equation  (21)  to  give 


^H.  Lamb.  Hydrodynamics.  Dover  Publications,  Inc..  Ne w  York  (1932). 
**/f  Goldstein.  Classical  Mechanics.  Addison-Wesley  Publishing  Co.. 
Reading.  MA  <19$0).  Ch  I 

1  D  Landau  and  E.  M.  Lifshitz.  Fluid  Mechanics,  Theoretical  Physics.  6, 
Addison  Wes  lev  Publishing  Co..  Reading.  MA  ( 1959),  5 4. 


The  decay,  when  present,  shifts  the  resonant 
frequency  so  that 


u>n  =  w2(n  =  0)  (22c) 

for  oj„  5:  0;  otherwise,  the  drop  does  not  os¬ 
cillate.  For  the  mode  n  =  2,  the  frequency 
o2  =  a)2/2n,  r2,  and  the  product  v2t2  were  calcu¬ 
lated  by  O’Konski  and  Thacher5  for  1-,  10-,  100-, 
and  1000- pm  water  droplets,  where  they  used 
y  =  72.0  dynes/cm  and  rj  =  0.884  X  10“2  poise. 

All  of  the  results  obtained  above  were 
previously  derived  by  various  techniques,  and  we 
have  shown  that  they  can  all  be  obtained  by 
methods  compatible  with  Rayleigh’s  original  ap¬ 
proach  to  the  problem.  Also,  as  we  have  shown, 
corrections  to  Rayleigh’s  results  are  more  trans¬ 
parent  by  this  approach. 

2.2  Electric  Field 

The  inclusion  of  an  electric  field  into  the 
analysis  presents  a  problem  of  considerable  com- 

‘  Lord  Rayleigh,  Proc.  Roy  Soc.  (London),  29  (18)9)  71;  The  Theory  of 
Sound.  II  MacMillan  Co.,  London  0896),  Ch  XX 

^H.  Lamb.  Hydrodynamics.  Doter  Publications,  Inc..  New  York  (1932). 

5C  T  OKonski  and  H.  Thacher.  Jr.,  J.  Phys.  Chem  .il  (19S3),  955. 
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plexity.  Therefore,  we  present  the  approach  and 
the  solution  to  the  problem  in  more  detail  than  in 
our  previous  discussion.  To  be  consistent  with  our 
preceding  analysis,  we  need  to  obtain  an  expres¬ 
sion  for  the  electromagnetic  energy  of  the  drop 
expressed  in  terms  of  the  a„(t)  in  equation  ( 1 ).  This 
energy  can  be  added  to  the  existing  Lagrangian, 
which  can  then  be  used  to  obtain  the  equations  of 
motion. 

A  convenient  form  for  the  electro¬ 
magnetic  energy  stored  in  the  drop  is  given  by 
Jackson:16 

UE  =  -1^-f  E‘  •  Edr  ,  (23) 

where  E'is  the  electric  field  inside  the  drop, 
E  is  the  field  in  the  absence  of  the  drop,  e  is  the 
dielectric  constant  of  the  drop,  and  the  integral 
covers  the  volume  of  the  drop.  Thus,  our  problem 
is  to  find  a  solution  for  E1  before  proceeding. 

For  a  dielectric  body,  we  have  from 
Maxwell's  equations 

V  X  E  =  0 

(24) 

V  •  D  =  e  V  •  E  =  0 


for  the  exterior  region  and 

=  Z  rnBnP„(cos  8)  (28) 

n 

for  the  interior  region  of  the  drop.  The  last  term  in 
equation  (27)  represents  the  potential  of  the  ex¬ 
ternal  electric  field,  E,  parallel  to  the  z-axis. 

The  applied  electric  field,  E,  can  be 

written 

E  =  E(f  cos  8  +  8  sin  8)  ,  (29) 

r  and  6  are  unit  vectors  in  their  respective  direc¬ 
tion,  and  from  equation  (25) 

Ef  =  -V^‘ 

Then  using  equation  (27),  we  have 
E  E1  =  -E  2  Bnrn-l[n  cos  8  Pn(cos  8) 

n  L 

-  sin  8  Pn(cos  0)J  .  (30) 

From  the  recursion  relation  for  Legendre  poly¬ 
nomials,  with  fx  =  cos  8,  we  have 

(!  -  M2)  =  nPn-l  -  n^n  , 


The  first  of  equation  (24)  implies  that  the  electric 
field  can  be  derived  from  a  scalar  potential,  or 

E  =  -v*  ,  (25) 

and  the  second  of  equation  (24)  implies  that 

=  0  .  (26) 

The  appropriate  solutions  to  equation  (26)  for  our 
problem  can  be  written 

oo  ^ 

*l'°  ~  Z  n-n  Pn(cos  8)  —  Er  cos  8  (27) 

n-0  r 


where  we  have  discontinued  writing  the  argument 
in  the  Legendre  polynomials.  Using  this  result  in 
equation  (30),  we  have 


E  E1  =  -E  Z  nBnrn_1Pn_, 


(31) 


and 


X2n  S'*  rx{6,  t)  _  _ 

d0Jb  sin^d 0jf,  r-’ E  •  E'dr 

=  ~2nE  rn+2Pn_,  dfi  . 


l)  Jack  ton.  Classical  Electrodynamics,  John  Wiley  and  Sons,  New  York 

fl975h  IbO  The  derivation  of  equation  (23)  is  not  trivial,  and  Jackson  yives 
considerable  attention  to  the  derivation 


Substituting  this  result  into  equation  (23)  we  have 

„  e  —  i  r  2b,  , ,  v  nB» 

Ue  4  E\~Y  3  +  n5.  n  +  2 


X  .  (32) 

where  we  have  separated  the  term  n  =1  from  the 
remainder  since  this  term  is  the  integral  over 
/x  of  r3,  which  is  proportional  to  the  volume  and  is 
the  constraining  equation  (3). 


To  evaluate  the  terms  in  the  sum  in 
equation  (32),  we  use  equation  ( 1 )  and  expand  the 
powers  of  ao  through  linear  terms  in  ak  or 

rn+2Pn_,  dM  -  /'  [ 

+  (n  +  2)aS+1  S'  akPkJ  P„- ,  d/x 


2(n  +2)  n+1 
(2n-l)a°  a"-' 


(33) 


for  n  >  1.  The  resulting  equation  (33)  can  be 
used  in  equation  (32)  to  give 


UE  = 


+ 


n  +  1 
2nTT 


Bn+ia"  2an 


’  (34) 


To  evaluate  Bn  and  An  in  equations 
(27)  and  (28),  we  need  to  apply  the  boundary 
conditions  on  the  fields  at  the  surface  of  the  drop. 
These  are  that  the  tangential  components  of  the 
electric  field  are  continuous,  E{  =  E°,  and  the 
normal  components  of  the  electric  displacement 
are  continuous,  D„  =  D„,  where  the  subscripts  t  and 
n  on  the  vector  components  denote  tangential  and 
normal,  respectively.  (The  use  of  n  here  and  in 
the  summing  index  need  cause  no  confusion  since 
no  sums  are  performed  over  the  field  components.) 
The  continuity  of  the  tangential  electric  field  can 
be  easily  shown  to  be  equivalent  to  the  continuity 
of  the  potential  so  that  i//‘  =  i j/°  at  the  surface. 
The_ boundary  condition  on  the  normal  component 
of  D,  however,  requires  some  consideration.  Since 
by  equation  (25)  E  =  —V  tfr  and  since  we  are 
assuming  that  D‘  =  t  E'  and  D°  =  E°  with 
n  a  unit  vector  normal  to  the  surface,  we  have 

e  n  •  v  =  n  •  v  l/'0  (35) 

at  the  boundary  of  the  drop.  To  construct  n,  we 
note  that  the  vector  dT  =  r  dr  +  Or  dd  lies  in 
the  surface  if  r  is  given  by  equation  ( 1 )  and  the 
unit  vector  0  in  the  0  direction  also  lies  in  the 
surface  since  we  are  assuming  that  the  drop  is 
symmetrical  about  the  z-axis.  A  vector  normal  to 
the  surface  can  be  formed  by  the  cross  product  of 
these  two  vectors  or 


and  we  have  shifted  the  summing  index  so  that  the 
lower  limit  on  n  is  1.  We  have  also  replaced  the 
product  ana£+2  by  anan+2  since  the  corrections 
to  ao  are  of  order  a£  so  that  corrections  to  this 
quantity  would  be  of  order  a£,  which  are  ne¬ 
glected.  Also,  in  equation  (33),  we  have  expanded 
r  only  through  terms  linear  in  ak  because  Bn  in 
equation  (34)  for  n  >  1  have  to  depend  on  ak 
since  they  vanish  for  undistorted  spherical  drops. 
(That  they  vanish  appears  in  the  subsequent  an¬ 
alysis.)  Therefore,  to  obtain  an  expression  valid 
through  terms  of  order  ajj,  it  is  necessary  to  obtain 
B,  through  terms  of  order  ajj  for  contributions 
from  the  first  term  in  equation  ( 34),  while  we  need 
only  first  order  terms  for  Bn  with  n  >  1.  To 
proceed,  we  have  to  return  to  the  problem  of 
evaluating  the  Bn  and  consequently  the  An  of 
equations  (27)  and  (28). 


diT  X  0 

!  dr  x  0! 


r  r  dd  -  0  dr 
(r2  d02  +  dr2)1'2 


(36) 


which  can  be  used  in  equation  (35)  to  obtain  the 
normal  derivative.  Thus,  to  reiterate,  the  boundary 
conditions  on  equations  (27)  and  (28)  are 


r  =  r 


d4>° 

dr 


(Jr\er 

\rd0/ r  dO 


(37) 


evaluated  at  the  boundary  with  r  given  by 
equation  ( 1 ). 
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To  keep  track  of  the  order  of  the  correc¬ 
tions,  it  is  convenient  to  rewrite  equation  (1)  as 


r  =  ao  +  S  2'  akpk 

and  expand  the  various  powers  of  r  in  powers 
of  5;  also,  we  assume  that 

A„  =  AS?  +  SK'  +  <$2Aln21  ...  , 


m2  ( 1 2mJ  +  2m2  —  1 8m  4-  10) 

(2m  +  l)2(2m  —  l)(2m  +  3)[e(m  —  l)  +  m] 


When  the  results  given  in  equation  (39)  are 
substituted  into  equation  (34),  we  have 


UE  =  - 


(e— l)E2a 


2(e  +  2) 


-■  +  ,, 


5(e  4-  2) 


Bn  ==  B‘n0’  +  «B5,'»  +  W  ...  .  (38) 

Then  in  the  final  result,  we  let  S  =  1  since  it 
serves  only  as  an  artifice  to  keep  the  terms  in 
order.  When  this  is  done  in  equation  (38), 
we  obtain 

A'n’  =  B„'  =  0  ,  n  *  1 


4~  V  -y  £  (fjm+2amam+2  +  Hma^,)  I 

E  +  2  m=l  J 

(40) 

To  form  the  full  Lagrangian  for  the  problem, 
we  need  only  add  UE  to  Us  so  that  L  = 
T  -  Us  —  UE,  and  the  equation  of  motion  for 
an  is  given  by  the  usual  procedure. 


A*0'  = 


(g  -  l)Ea^ 
e  4-  2 


^  <9L  _  .  (41) 

dt  <?an  dan  da„ 


B<°>  =  - 


3E 
£  4-  2 


The  result  of  doing  this  is  to  add  to  the  right  side  of 
equation  (22)  the  term  —  dUE/dan  so  that 


3(e  —  1)E  n(n -f  lXe  -  l)an+, 
~l+2 “  (2n  +  3)[(e+  l)n+  lj 


+  -^i.  35+' 
2n-  1  * 


B'n"  = 

and 


B«2' 


3(e  -  1)E 

(n  +  lX2n  4-  l)an+. 

e  +  2 

(2n  +  3)[(f+  l)n+  lJaS 

9(r  —  1)E 

(e  +  2)2aii 


47rpa' 
n(2n  +  1) 


an  +  8nrj 


a(n  —  1)  . 

ai 


4Try(n  -  l)(n  +  2) 

+  - 2rT  4-~l - 3" 

3(e  -l)2E:a  x  3(c  —  l)2E2a 

5(r  4-  2)2  "•*  (e  +  2)2 

X  (Cn+2an+2  Gnan_2  4-  2Hnan)>  (42) 

which  are  the  equations  of  motion  for  a  drop  in  an 
electric  field.  The  right  hand  side  of  equation  (42) 
appears  to  be  new. 


X  X(Gmamam-2  +  Hmaj^  i  (39) 
where 

m(m  —  l)(2m  —  l)(e  +  2) 

°m  ~  (2m  —  3)(2m  +  l)l£(m-  l)  +  m]  ’ 

m[r(m  -  1)(  12m3  4-  10m2  -  12m  -  l)] 
(2m  4-  l)2(2m  -  lX2m  +  3)[e(m  -  1)  +  m] 


Several  aspects  of  equation  (42)  are  as 
should  be  expected  from  the  general  symmetry  of 
the  problem.  That  is,  the  right  hand  side  is  depen¬ 
dent  on  the  square  of  the  electric  field  so  that  a 
reversal  of  the  field  does  not  alter  the  results. 
Further,  the  n  =  2  mode  is  the  only  mode  driven 
directly  by  the  electric  field,  and  this  couples  to 
the  n  =  4  mode  and  consequently  to  higher 
modes  with  n  even.  Consequently  if  no  further 
perturbation  couples  the  odd  and  even  modes,  only 
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those  modes  for  n  even  need  to  be  considered. 
Some  of  these  results  would  have  been  immedi¬ 
ately  obvious  if  we  had  used  Maxwell’s  stress 
tensor  to  evaluate  the  forces  on  the  drop,  but  then 
the  normal  mode  result  of  equation  (42)  would 
have  been  lost.  If  we  ignore  the  coupling  to  the 
higher  modes  in  equation  (42),  the  resonant  fre¬ 
quency  of  the  nth  mode  is  given  by 


(ti 


2 

n 


n(n-  1) 
pa3 


X 


yin +  2) 


6(e  —  l)2E2a(2n  4-  l)Hn  ~| 
4jT(n-  l)(e  +  2)2  J  ■ 

(43) 


This  result  is  similar  in  appearance  to  the  result 
derived  by  Rayleigh  for  a  charged  drop, 


iO 


2 

n 


n(n-  1) 
pa3 


y(n  +  2) 


3.  DISCUSSION 

The  result  given  in  equation  (42)  is  to  this  point 
general  in  that  no  assumptions  have  been  made 
concerning  the  nature  of  the  electric  field  other 
than  those  implied  in  equations  (24)  and  (25).  In 
this  section,  we  discuss  results  from  equation  (42) 
for  three  different  types  of  fields:  a  static  electric 
field,  an  alternating  electric  field,  and  an  amplitude 
modulated  high  frequency  field. 

3. 1  Static  Electric  Field 

In  the  case  of  a  charged  drop  under  the 
influence  of  its  self-electric  field,  the  modes  of 
oscillation  are  uncoupled  so  that  a  simple  criterion, 
wn  =  0,  is  usually  taken  as  the  onset  of  insta¬ 
bility.  In  the  present  case  given  in  equation  (42), 
we  see  that  the  modes  are  coupled  so  that  the 
condition  for  instability  created  by  a  large  electric 
field  becomes  complicated.  If  we  assume  that 
an  and  an  in  equation  ( 42)  vanish,  then  we  can  write 


with  the  charge  on  the  drop,  Q  =  Ea2.  In  Ray¬ 
leigh’s  result,  the  frequency  decreases  as  the 
field  (charge)  increases,  but  in  equation  (43)  this 
condition  occurs  only  if  H„  >  0.  The  condition 
for  no  shift  in  the  frequency  in  equation  (43)  is 
Hn  =  0,  which  for  the  first  few  modes  are 
n  =  2,  e  =  2.846;  n  =  4,  e  =  1.588;  n  =  6, 
e  =  1.346;  and  n  =  8,  e  =  1.  Perhaps  a  fair 
comparison  with  Rayleigh’s  result  is  for  a  con¬ 
ducting  drop,  which  can  be  obtained  from  equation 
(43)  by  letting  the  dielectric  constant  be  infinite  or 


“>i i(f  =  °°)  - 

n(n  -  1) 
pa' 

|y(n  +  2)  - 

6E2a 

4n 

n(  1 2n’  +  1  On2  —  1 2n  -  1 )  1 
(n  -lX2n  +  l)(2n  -l)(2n  +  3)J  '  ' 

This  resulting  equation  (44)  shows  the  same 
behavior  as  the  result  of  Rayleigh,  but  with  a  more 
complicated  dependence  on  mode  number,  n. 


Dnxn  j  ^y^n.2 

+  3Xy(Gn+2xn+J  +  Gnxn_2)  ,  (45) 

where 


Dn  =  An  -  6\yHn 

4nfn  —  1  )fn  +  2) 


with  Gn  and  Hn  given  in  equation  (39).  For  low 
fields  (y  ~  0),  the  coupling  between  modes  can 
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generally  be  ignored;  but  for  larger  fields,  which 
would  be  required  for  the  drop  to  become  unstable, 
the  coupling  cannot  be  ignored.  Thus  to  investi¬ 
gate  the  instability  in  the  system  of  equations  in 
equation  (45),  it  is  necessary  to  consider  the 
coupled  equations  in  detail. 

If  we  write  equation  (45)  in  matrix  nota¬ 
tion,  then 


Tx  =  F  ,  (46) 


where 

Tnn  =  Dn5n'n  —  3XGn5n.n+2  —  3AGn5n'n_ 2 

and  x  and  F  are  column  matrices  (vectors).  The 
matrix  F  consists  of  the  single  element,  (3/5)Ay. 
The  formal  solution  to  equation  (46)  is  given  by 

x  =  T‘F  ,  (47) 

where  T“'  is  the  inverse  of  the  matrix  T.  The 
inverse  of  the  matrix  T  contains  the  determinant 
of  T  in  the  denominator,  and  when  the  deter¬ 
minant  vanishes,  the  solution  equation  (47)  be¬ 
comes  unstable.  Thus  the  lowest  value  of  electric 
field  (smallest  positive  y)  for  which  the  deter¬ 
minant  vanishes  gives  the  onset  of  instability  of  the 
drop.  Since  the  matrix  T  has  single  off-diagonal 
elements,  it  is  a  simple  matter  to  obtain  a  recursion 
relation  for  the  determinant,  A.  If  the  matrix  is 
truncated  to  contain  N  terms  (Tnn’  —  0, 
n  or  n'  >  2N),  then 

A2n  =  D2NA2n-2  -  (3XyG2N)2A2N_4  -  (48) 

with  A0  =  1  and  A2n-2  =  0  for  N  negative. 
The  result  given  in  equation  (48)  can  be  used  to 
obtain  the  smallest  positive  value  of  y  such  that 
A,N(y)  =  0.  The  first  two  of  these  are  from 
equation  (48), 


A2(y)  =  A2  -  6AyH2  ,  A2(y)  =  0  gives 


y2  = 


a2 

6AH, 


(49) 


and 


A4(y)  =  (A4  -  6AyH4) 

X  (A2  -  6AyH2)  -  (3AyG4)2 


’  (50) 


H4A2  +  H2A4 
y 4  3A(4H4H2  -  G4) 

[(h4a2  -  H2A4)2  +  a2a4g2]12 
3A(4H4H2  -  G2) 


which  are  the  simplest  roots  to  obtain  algebraic 
expressions  for.  The  ambiguity  of  sign  in  equation 
(50)  was  determined  so  that,  y4  given,  there  was 
the  smallest  positive  root  in  the  range  of  dielectric 
constants,  1  <  e  <  °°.  The  result  for  y2  in 
equation  (49)  is  not  physical  (y2  <  0)  for  e  < 
52/37  and  becomes  infinite  at  t  =  52/37  (H2 
vanishes)  and  is,  in  general,  not  even  a  good  first 
estimate.  This  result  is  not  surprising  since  y2 
totally  ignores  the  coupling  between  modes.  The 
result  for  y4  given  in  equation  (50)  is,  on  the  other 
hand,  a  good  estimate  even  for  small  values  of  e. 
The  lowest  values  of  y  for  which  A2N(y)  =  0 
have  been  calculated  for  a  few  values  of  N  and  an 
extended  range  of  e  and  are  given  in  table  1 .  The 
value  of  y4  differs  from  the  value  of  y,„  only  by 
16percentfor  e  =  1.1;  for  water  (e  =  78.2),  y4 
differs  from  yl0  by  1.3  percent,  and  for  larger  e 
the  difference  is  approximately  1  percent. 


TABLE  1.  CRITICAL  VALUES  OF  y  (E2a/y) 
ARE  LOWEST  VALUES  OF  y  FOR  WHICH 
Ak(yk)  =  0  FOR  k  =  4,  6,  8,  AND  10 


e 

y  6 

y« 

yio 

i.i 

7009 

6109 

6054 

6049 

1.3 

817.6 

728.2 

720.5 

720.0 

1.5 

310.9 

280.4 

278.2 

278.1 

2.0 

90.16 

83.71 

82.80 

82.78 

5.0 

12.87 

12.42 

12.40 

12.39 

10.0 

6.511 

6.367 

6.362 

6.362 

50.0 

3.603 

3.556 

3.555 

3.555 

78.2 

3.406 

3.364 

3.363 

3.363 

100 

3.332 

3.302 

3.291 

3.291 

oo 

3.078 

3.044 

3.044 

3.044 
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Upon  comparing  the  results  in  table  1 
with  the  results  obtained  by  Brazier-Smith  et  al,17 
we  find  a  discrepancy.  The  value  of  E(a/y)l/i 
=  l  .625  (y  =  E2a/y)  given  by  them  (along  with 
references  to  previous  work)  corresponds  to  our 
value  of  1.745  for  e  =  °° .  For  water,  which  is 
used  in  their  experiment,  we  obtain  1.834.  This 
discrepancy  is  caused  possibly  by  the  assumption 
that  the  drop  is  constrained  to  an  ellipsoid  of 
revolution  in  their  treatment  of  the  problem. 

The  amplitude  of  the  x2  mode  is  easily 
obtained  by  using  equation  (45)  along  with  the 
determinant  in  equation  (48)  to  give  for  the  2Nth 
approximate  solution 

P;M 

x2(2N)  =  ,  (51) 

a2n 

where  the  P2n  obey  the  same  recursion  relation 
as  AjN  given  in  equation  (48).  The  initial  values 
of  PjN  are 

P2  =  3Xy/5  (52) 


Figure  1.  Amplitude  of  second  mode,  X2,  as  function 
of  y  (E2a/y)  for  dielectric  constant 
f  =  78.2  (A)  and  f  =  -  (B). 


and 


P4  =  3XyD4/5  ,  (53) 


which  are  sufficient  to  generate  all  Pjn  from 
equation  (48).  Having  determined  the  value  of 
x2(2N),  we  can  generate  all  the  x„(2N)  for 
n  5  2N  by  using  equation  (45)  and  the  condition 
x„  =  0.  The  variation  of  x2(lO)  as  a  function  of  y 
is  shown  in  figures  I  and  2  for  two  values  of 
dielectric  constant,  e  =  78.2  and  e  =  °° .  The 
amplitudes  of  the  xn  mode  are  approximately  an 
order  of  magnitude  larger  than  the  xn+2  mode  in 
the  moderate  y  region.  All  the  modes,  x„,  di¬ 
verge  at  the  smallest  value  of  y,  satisfying  the 
equation  A,0(y)  =  0.  When  the  xn  become 
large,  the  results  are  questionable,  and  a  possible 
measure  of  the  range  of  xn  that  are  reasonable  is 
from  equation  (5)  (with  xn  =  a„/a). 


«0 

a 


=  1 


2n  +  1 


(54) 


1 ' P  R.  Brazier Smith.  M  Brook.  J.  Lantham.  C.PR ■  Sounder*,  and  M.  H 
Smith.  Pfytt-  Hoy  Sot-  (l.omjon).  A 322  (1^7/).  $23. 


Figure  2.  Amplitude  of  fourth  mode,  xi,  as  function 
of  y  (EJa/y)  for  dielectric  constant 
f  =  78.2  (A)  and  f  =  -  <BI). 
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When  ao/a  becomes  small,  the  entire  results  be¬ 
come  questionable  because  all  the  expansions 
used  in  the  derivation  equation  (45)  assume 
ao  >  (  ak  i  implicitly. 


If  the  drop  were  distorted  in  the  shape  of 
an  ellipsoid  of  revolution  along  the  direction  of  the 
field,  we  could  write 


r2  sin2  6  r2  cos2  6  _  ,  iecs 
A2  B2  '  (55) 

If  we  let  6  —  tt/2,  we  have  from  equation  ( 1 ) 
A,  =  r(~)  =  £  anPn(0)  (56) 


or 


A  =  x2n(-l)"(2nn)  yj  .  (57) 


where  A  =  A,/a 

and^n^  is  a  binomial  coefficient. 

If  we  let  8  —  0  or  n  in  equation  (55),  then  from 
equation  ( 1 ) 


B,  = 


dfi  =  0)  =  ao  +  £'  an  , 
r[8  =  jr)  =  a„  +  £'  (—  l)nan 


Then 


B  -  y  +  *2"  *  <58) 

where  B  =  B,/a.  If  the  distortion  were  an  ellip¬ 
soid  of  revolution  as  given  in  equation  (55),  the 
volume  of  the  drop  would  be  given  by  A2B  =  1 . 
For  e  =  ~  .  ao/a  was  computed  from  equation 
(54),  A  by  equation  (57),  and  B  by  equation 
(58)  as  a  function  of  y,  and  the  results  are  shown 
in  Figure  3.  A  measure  of  how  closely  the  drop 
approximates  an  ellipsoid  of  revolution  is 
1  -  A2B;  for  small  values  of  y,  this  quantity  is 
quite  small.  However,  for  y  =  2.8,  1  —  A2B~ 
0.28.  and  ao/a  is  still  large  (~0.9).  The  ellip¬ 


soidal  approximation,  when  viewed  in  this  light, 
would  not  seem  very  good  and  perhaps  would 
indicate  the  source  of  the  difference  in  the  critical 
field  values  obtained  here  and  elsewhere  by  earlier 
workers.  (See  fig.  3  on  p.  16.) 

3.2  Dynamical  Equation  for  Varying 
Electric  Fields 

The  dynamical  equation  (42)  is,  to  our 
knowledge,  a  new  result  in  the  sense  that  the  losses 
due  to  viscosity  and  a  finite  dielectric  constant 
have  been  included  in  the  analysis.  Further,  the 
entire  equation  of  motion  was  derived  in  a  con¬ 
sistent  manner  by  using  Rayleigh’s  original 
method.  For  convenience,  we  can  write  equation 
(42)  in  the  form 

3 

Mnxn  +  Mnrnx„  +  Mnaif,xn  =  ,  AE2atfn  2 

,  ,  <59> 

+  3XE  a(Gn+2xn+2  +  Gnxn_2)  , 

where 


Mn  = 


47rpa' 


n(2n  +  1) 

(M,  is  the  total  mass  of  the  drop), 
2rja(n  -  l)(2n  +  1) 


rn  = 


pa 


urf  =  u)2(0)-  n(2n  +  1)H„ 

where  u>n(0)  is  the  frequency  in  the  absence  of  an 
electric  field  given  in  equation  (22a)  and  H„  is 
given  in  equation  (39).  The  quantity  X  is  defined 
in  equation  (45).  If  the  electric  Field  is  assumed  in 
the  form 

E  =  E0  cos  u>t  ,  (60) 

then 

E2  =  -y  +  -y  cos  2wt  (61) 
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Figure  3.  Quantities  %/a,  a2/a,  and  1  —  A2B  obtained  from  normal  mode  analysis  versus  field  strength  y. 


Thus,  the  driving  force  in  equation  (59)  for  the  and 
n  =  2  mode  varies  at  twice  the  frequency  of  the 
electric  field,  and  resonance  for  small  amplitude  of 
oscillation  should  occur  near  2to  =  wn.  ^  _ 

If  we  ignore  the  coupling  to  higher  modes  where 
(Xo  =  0)  in  equation  (59)  and  replace  E2  by 
El/ 2  in  we  have  R  = 


3AE2a 

10M,R, 


cos  (2(ut  -  02) 


R2  =  [(tof  -  4w2)2  +  4u>2fl]' 


x2  +  r2x2  +  w2x2 


3\Ena  JAEpa  -  ,,  o>T2 

“  TomT  +  TomT  “*  2“'  •  m  *■  =  4  -  w 

where  we  have  used  the  results  of  equation  (61 ).  Also,  in  the  derivation  of  equation  (63),  we  have 

The  result  given  in  equation  (62)  is  a  standard  ignored  terms  of  higher  power  than  Ep  in  the 

linear  equation,  and  the  solution  can  be  written  electric  field. 
x2(t)  =  x°2  +  <2(t)  with 

2  Thus  the  results  given  in  equation  (63) 

o  _  3AE0a  show  that  in  small  electric  fields  the  drop,  in  the 


10M2to2 


fundamental  mode,  is  similar  to  an  ordinary  reso- 
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nant  circuit  with  damping;  however,  the  driving 
field  on  the  drop  has  a  frequency  double  that  of  the 
applied  field.  The  higher  modes  are  driven  by  the 
electric  field  indirectly  by  coupling  through  the 
lower  modes  as  shown  in  equation  (59),  and 
approximate  solutions  can  be  obtained  by  a  pertur¬ 
bation  theory  solution  of  the  equations  of  motion. 
If  the  coupling  is  ignored,  the  amplitude  of  vibra¬ 
tion  of  the  nth  mode  can  be  written 

x„(t)  =  atf'e  r"/2  cos  (flnt  +  an)  (64) 
for 

fS/4  <  u>n 


and 


with  ton  and  r„  given  by  equation  (59).  The  re¬ 
sult  in  equation  (64)  shows  that  the  nth  mode 
decays  more  rapidly  than  the  fundamental  (n  =  2); 
that  is,  since 

(n  ~ I )(2n  +  I)  „ 

>n - ;  l2 


that  in  certain  cases  the  resonant  frequency  could 
be  unshifted  if  the  viscosity  were  small  enough. 
Care  must  be  exercised  in  extending  equation  (66) 
to  higher  electric  fields  because  the  coupling  to 
higher  modes  given  in  equation  (59)  can  have  more 
important  consequences.  Also,  for  higher  electric 
field  strengths,  the  term  in  w„  in  equation  (59) 
involving  the  electric  field  has  a  time  dependent 
part  that  could  have  some  interesting  conse¬ 
quences  in  large  electric  fields. 

The  dynamical  equation  (62)  or  the  more 
general  equation  (59)  can  be  used  in  more  compli¬ 
cated  situations  where  the  electric  field  is  of  an 
impulse  nature  such  as  that  caused  by  the  presence 
of  charged  drops.  A  number  of  sources  of  electrical 
and  mechanical  disturbances  have  been  consid¬ 
ered  by  Brook  and  Lantham  in  their  study  of 
modulation  of  radar  echo  from  rainstorms.18 

A  second  and  perhaps  a  more  interesting  case 
of  varying  fields  is  the  electric  field  as  an  amplitude 
modulated  high  frequency  wave.  The  carrier  may 
be  radar  or  an  infrared  laser,  and  the  amplitude 
modulation  can  be  chosen  near  a  resonance  of  the 
drop.  Such  a  field  can  be  represented  by 

E  =  E0(l  +  mo  cos  wt)  cos  w0t  ,  (67) 


the  n  =  4  mode  decays  approximately  six  times 
faster  than  the  fundamental.  The  shift  in  frequency 
due  to  viscosity  and  electric  field  can  be  written 

Awn  =  w„(0)  -  On  .  (65) 

Assuming  the  electric  field  and  the  viscosity  small, 
this  can  be  written 

Awn 

w„(0) 

-jsSoi  +  mSt 


where  w0  >>  w,  and  particular  space  depen¬ 
dence  has  been  ignored  because  we  assume  that 
the  drop  radius  is  small  compared  with  the  wave¬ 
length,  \>,  of  the  carrier  (\,  =  2ttc/w0).  We 
need  the  square  of  the  field,  which  is  given  by 

E2  =  ^  ^1  +  ^  +  2 mo  cos  wt 

+  cos  2wt  j  (68) 


(66) 

In  general,  the  shift  in  frequency  given  in  equation 
(66)  is  positive.  However,  Hn  given  in  equation 
( 39)  can  be  positive  or  negative  depending  on  E  so 


If  we  ignore  all  the  high  frequency  terms  varying  at 
frequencies  near  2w0  and  if  the  depth  of  modula¬ 
tion  mo  is  not  large,  the  terms  involving  ms  can 


Brook  and  D.  J.  Lanlhan,  J.  Geophys.  Res.,  7.3  ( 1968).  7137. 
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be  ignored.  Substituting  equation  (68)  into  equa¬ 
tion  (59),  we  obtain 

.  ,  -  •  ,  ,  3AEoa 

xn  rnXn  +  wnXn  —  — 

X  (1  +  2mo  cos  wt )Sn2+  — 

X(l  +  2mo  cos  tat)(Gn+2xn+j  +  Gnxn_2)  , 

(69) 

where  as  in  equation  (62)  we  assume  that  the 
term  E2  in  is  replaced  by  Eo/2.  Also  im¬ 
plicit  in  the  equation  is  the  assumption  that  the 
dielectric  constant,  e,  wherever  it  appears,  is  to 
be  evaluated  at  the  carrier  frequency,  oi0.  Only  in 
cases  of  extreme  dispersion  is  this  invalid,  and 
such  cases  have  to  be  treated  quite  differently.  The 
solutions  to  equation  (69)  are  the  same  as  in 


equation  (62)  for  n  =  2,  except  that  the  fre¬ 
quency  here  is  at  half  the  value  used  in  equation 
(62). 

The  time  varying  displacement  of  the  drops 
could  be  studied19  by  amplitude  modulating  a  high 
powered  laser  and  observing  the  modulation  of  the 
reflected  signal  by  a  second  low  power  laser.  Some 
of  the  analysis  given  by  Brook  and  Lantham18 
would  apply  to  this  case,  except  that  the  modula¬ 
tion  of  the  reflected  signal  from  the  low  energy 
laser  would  be  at  the  frequency  u>  of  equation 
(60),  rather  than  a  distribution  of  different  fre¬ 
quencies.  By  observing  the  reflected  signal  of  the 
low  energy  laser  and  sweeping  through  a  range  of 
values  of  to,  one  could  obtain  some  estimate  of 
the  distribution  of  particle  size. 

Brook  and  D.  J.  Lantham.  J.  Geophys.  Res..  73  (1968).  7137. 

E.  Wonman,  Possible  Use  of  Two  Laser  Beams  to  Determine  Panicle- 
Sue  Distribution,  Harry  Diamond  Laboratories  HDL-TR- 18 78 (January  1979). 
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